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Abstract 

Unsteady transonic flow calculations for 
aerodynamical ly Interfering airfoil configura- 
tions are performed as a first-step toward 
solving the three-dimensional canard-wing inter- 
action problem* These calculations are 
performed by extendi ng the XTRAN2L two- 
dimenslonal unsteady transonic small -disturbance 
code to include an additional airfoil. Unsteady 
transonic forces due to plunge and pitch motions 
of a two-dimensional canard and wing are 
presented* Results for a variety of canard-wing 
separation distances reveal the effects of aero- 
dynamic Interference on unsteady transonic 
airloads. Aeroelastlc analyses employing these 
unsteady airloads demonstrate the effects of 
aerodynamic interference on aeroelastlc 
stability and flutter. For the configurations 
studied, increases In wing flutter speed result 
with the Inclusion of the aerodynamlcally inter- 
fering canard. 


Nomenclature 


a 

b 


c 



D 


h 


h l 

k 


m 

M 


q 

<lF 

Qij 


r a 

s 

S 


t 

u 


Uf 

Vf 


nondimensional elastic axis location, 
positive aft of mid-chord 
wing semi -chord 
wing chord 

canard lift coefficient 

wing lift coefficient 

steady pressure coefficient 

critical pressure coefficient for sonic 

flow 

canard-wing vertical separation 
distance in units of wing chord (see 
Fig. 1(a)) 

plunge displacement in units of wing 
chord 

plunge pulse amplitude 
wb/U, reduced frequency 
mass of the wing per unit span 
freestream Mach number 
(U/(bw a /IT)) 2 . nondimensional flight 
dynamic pressure 
(Up/(bco a /T7) ) 2 , nondimensional 
flutter dynamic pressure 
unsteady aerodynamic coefficient, 
ith generalized force due to the 
jth mode of motion (see Table 1) 
radius of gyration of wing about 
elastic axis 

o+l w, Laplace transform variable 
canard-wing horizontal separation 
distance in units of wing chord (see 
Fig. 1(a)) 
time 

freestream velocity 
flutter speed 
U F /(bw a /ir), flutter speed 
Index 


x distance from leading edge positive 

downstream along canard or wing 
x a distance from elastic axis to mass 

center nondimensionallzed by b 
{z} state vector 

a angle of attack of canard or wing 

wing pitch angle 

o 0 mean angle of attack of canard or 

wing 

pitch pulse amplitude 
ACp lifting pressure coefficient 

AM Mach number contour Interval 

p m/npb 2 , wing mass ratio 

p freestream air density 

t Ut/b, nondimensional time 

w angular frequency 

u)h uncoupled wing plunge natural 

frequency 

<u a uncoupled wing pitch natural 

frequency 


Introduction 

Computational methods employing linear 
theory have been developed for the prediction of 
unsteady flowflelds about aerodynamlcally Inter- 
fering lifting surfaces. 1 * 2 These methods are 
extensions of linearized theory for single lift- 
ing surfaces to more complicated configurations 
such as two lifting surfaces In tandem. Appli- 
cations of these methods are restricted to 
subsonic or supersonic flows, though, because of 
the underlying linear theories on which they are 
based. In the transonic regime, computational 
methods for modeling aerodynamic Interference 
flowfields about oscillating multiple lifting 
surface configurations have yet to be developed. 

Steady transonic flowflelds and Inter- 
ference effects about two-dimensional canard- 
wing systems have been studied by Shankar, 
Malmuth, and Cole. 3 Using transonic small-dis- 
turbance theory, the Interference flowflelds 
were computed in preparation for solving the 
three-dimensional interaction problem. In Ref. 
3, a double grid arrangement was used whereby 
the canard and wing were placed In separate 
computational domains. Results showed a favor- 
able increase in overall lift established by 
appropriate placement of canard and wing. 
Steady transonic computational results for 
three-dimensional canard-wing configurations 
Including comparisons with experimental data 
have been reported in Refs. 4 and 5. Shankar 
and Malmuth 1 * presented computations for a few 
canard-wing configurations obtained by placing 
the two surfaces in a sheared fine grid system 
that Is embedded in a global Cartesian crude 
grid. A weakening of the wing shock due to the 
canard downwash was illustrated by Mach number 
contour plots. Shankar and GoebeV 5 developed a 
local numerical mapping procedure where the 
leading and trailing edges of the two surfaces 
are treated as constant coordinate lines In the 
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computational domain. Computed results for a 
closely-coupled canard-wing research model were 
in good agreement with experimental data. Also 
in Ref. 5, inaccuracies were demonstrated in 
results computed using double grid arrangements 
for closely-coupled canard-wing systems. 

The purpose of this paper is to present 
unsteady transonic results for two-dimensional 
canard-wing configurations as a first-step 
toward solving the unsteady three-dimensional 
interaction problem. The objectives of the 
study were to: (1) develop an unsteady transonic 
computational capability for aerodynamical ly 
interfering airfoils; (2) investigate the 
effects of two-dimensional canard-wing aero- 
dynamic interference on transonic steady pres- 
sures and unsteady forces; and (3) determine the 
effects of canard-wing separation distance on 
transonic unsteady airloads, aeroelastic 
stability, and flutter, for a limited number of 
configurations. In this study, the aerodynamic 
calculations were performed by extending the 
XTRAN2L 6 transonic small -disturbance code to 
allow the treatment of an additional airfoil. 
The leading airfoil has been termed the "canard" 
and the trailing airfoil has been termed the 
"wing". Steady transonic pressure distributions 
and Mach number contour plots are presented for 
isolated airfoil and canard-wing configurations. 
Unsteady transonic forces due to plunge and 
pitch motions of the canard and wing are shown. 
Results for a variety of canard-wing separation 
distances were obtained to determine the effects 
of aerodynamic interference on unsteady tran- 
sonic airloads. These transonic airloads are 
used in aeroelastic analyses to demonstrate 
application of the XTRAN2L canard-wing computa- 
tional capability and to investigate the effects 
of aerodynamic interference on aeroelastic 
stability and flutter. 

Computational Procedures 
Transonic Code XTRAN2L 

The original LTRAN2 code of Ballhaus and 
Goorjian 7 was developed to time-accurately 
integrate the low-frequency transonic small - 
disturbance (TSD) equation with steady-state 
airfoil and wake boundary conditions. Houwink 
and Van der Vooren 8 extended the range of use- 
able frequencies by adding the time-derivative 
terms to the airfoil and wake boundary condi- 
tions. The resulting code was termed LTRAN2- 
NLR. The XTRAN2L code is an extensive modifica- 
tion of LTRAN2-NLR which was developed at NASA 
Langley Research Center. The program solves the 
complete TSD equation and includes monotone 
differencing, nonreflecting farfield boundary 
conditions, an improved grid, a pulse transient 
capability, and time-marching aeroelastic 
analyses. Details of the XTRAN2L algorithm 
development and modifications are given by 
Whitlow. 6 Details of the grid development and 
pulse capability are given by Seidel, Bennett, 
and Whitlow. 9 

Two-dimensional canard-wing calculations 
were performed by extending the alternating- 
direction implicit (ADI) solution procedure of 
XTRAN2L to include an additional airfoil. The 
program is now capable of computing unsteady 
transonic flowfields about relatively general 


interfering airfoil configurations. The present 
program coding, though, does not allow for over- 
lapping or co-planar configurations. 

The two-dimensional canard-wing geometry 
considered is shown in Fig. 1(a). Horizontal 
separation distance is defined in units of wing 
chord from wing midchord to canard midchord by 
S; vertical separation distance is defined In 
units of wing chord by D. The XTRAN2L grid near 
the canard-wing system Is shown in Fig. 1(b). 
Flat horizontal wakes are assumed for both 
canard and wing. The grid is based on the 
XTRAN2L default grid described in Ref. 9 and, 
for the configuration shown (S = 1.0, D s 0.25), 
has 116 and 66 points In the horizontal and 
vertical directions, respectively. In the 
streamwise direction, thirty and fifty equl- 
distantly spaced points are distributed along 
the canard and wing, respectively. One addi- 
tional point is added near the leading edge of 
each airfoil for better resolution. Grid 
stretching to the farfield is identical to that 
described in Ref. 9. For example, upstream of 
the canard and downstream of the wing, grid- 
points are stretched twenty wing chordlengths to 
the farfield boundaries. Above the canard and 
below the wing, gridpoints are stretched twenty- 
five wing chordlengths to the farfield bound- 
aries. Between the two airfoils, vertical grid 
lines are distributed equidistantly between the 
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(a) two-dimensional canard-wing 
geometry. 



(b) XTRAN2L grid near canard-wing 
system. 

Fig. 1 Geometry definition and finite- 
difference grid; 
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first gridpoint downstream of the canard trail- 
ing edge and the wing leading edge. Horizontal 
grid lines are distributed equidistantly between 
the airfoil centerlines with one additional line 
placed symmetrically below the canard centerline 
and another additional line placed symmetrically 
above the wing centerline. For different 
canard-wing separation distances, points are 
added or removed to maintain similar grid 
spacing. 

Pulse Transfer-Function Analysis 

Unsteady aerodynamic forces are computed 
for two modes of airfoil motion: pitch about the 
quarter chord and vertical translation (plunge). 
Typically, unsteady aerodynamic forces are 
determined by calculating several cycles of 
forced harmonic oscillation with the last cycle 
providing the estimate of the forces. Alterna- 
tively, harmonic forces may be obtained 
indirectly from the response due to a smoothly 
varying exponentially shaped pulse. 9 In this 
procedure, the airfoil is given a small 

prescribed pulse in a given mode of motion 
(either plunge or pitch) and the aerodynamic 
transients calculated. For pitch motion, the 
input pulse is given by 

. -0.25(t-17.5&t) 2 , n 

a = olq + a\e v ' (1) 

and for plunge motion, the input pulse is given 
by 

h ■ hle -0.25(,-17.5 a ,) 2 (2) 

where At is the nondimensional time step. The 
harmonic response is obtained by a transfer- 
function analysis using fast Fourier transforms 
(FFT). Use of the pulse transfer-function 
technique gives considerable detail in the 
frequency domain with a significant reduction in 
cost over the alternative method of calculating 
multiple oscillatory responses. 

Pulse transient calculations were performed 
using 1024 time steps with At set equal to 
5n/32. Plunge pulse and pitch pulse amplitudes 
were hi = 0.01 and 04 = 0.1°, respectively, 
for both canard and wing. 


Aeroelastic Model 

For aeroelastic analysis, the wing was 
assumed to have plunge and pitch degrees-of- 
freedom and the canard was assumed to be motion- 
less. Thus the canard influences the stability 
of the wing through aerodynamic coupling only. 
The structural equations of motion used for the 
wing are the classical equations for an airfoil 
section oscillating with plunge and pitch 
motions. 10 Aeroelastic parameter values select- 
ed were a = -0.5, x a - 0.2, r a = 0.5, 
wh/uiot = 0.3, and y = 60. 

Aeroelastic stability calculations were 
performed using a state-space aeroelastic model 
termed the PadG model. 11 * 12 This model is 
formulated by curve-fitting the unsteady aero- 
dynamic forces by a Padd approximating func- 
tion. 11 These approximating functions are then 
expressed as linear differential equations 
which, when coupled to the structural equations 


of motion, lead to the first order matrix 
equation 


U} s [A]{Z} (3) 

where (Z) is a state vector containing displace- 
ments, velocities, and augmented states. The 
augmented states model the unsteady airloads. 
Stability analyses are performed by a linear 
eigenvalue solution of the Padd model. Real and 
imaginary parts of the eigenvalues (damping and 
frequency, respectively) are plotted in a 
dynamic pressure "root- locus” type format In the 
complex s-plane. 


Results and Discussion 


Steady pressure calculations for a simple 
case of two flat plate airfoils were performed 
first, to assess the XTRAN2L code modifications 
by comparison with an independent vortex-lattice 
method (VLM) program. The freestream Mach 
number was M = 0.5; the mean angle of attack was 
a 0 = 1.0° for both airfoils; the canard chord- 
length was selected to be 60% that of the wing 
chordlength; the horizontal and vertical separa- 
tion distances were S - 1.0 and 0 = 0.25, 
respectively. The lifting pressure coefficient, 
ACp, is plotted in Fig. 2 for isolated air- 
foils and closely-coupled canard-wing configura- 
tions. The XTRAN2L results are in excellent 
agreement with results from the VLM program. As 
shown in Fig. 2, the area between the Isolated 
airfoil AC p curves and the ( canard-wing AC p 
curves represents the steady aerodynamic inter- 
ference between the two airfoils. For the 
configuration shown, the canard induces a down- 
wash on the wing thus decreasing its AC p and 
lift. Conversely, the wing Induces an upwash on 
the canard which increases its AC p and lift. 


XTRAN2L / ISOLATED AIRFOILS 

XTRAN2L / CANARD-WING 

o VLM / ISOLATED AIRFOILS 



Fig. 2 Comparison of lifting pressure 

coefficients computed using XTRAN2L 
and a vortex lattice method (VLM) 
for flat plates at M = 0.5 and 
a 0 = 1.0° (S - 1.0, D = 0.25). 
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For transonic XTRAN2L computations, the 
NACA 0010 13 airfoil section was selected for 
both canard and wing. The freestream Mach 
number was set equal to 0,76; the mean angle of 
attack was a 0 = 1.0° for both airfoils; the 
canard chordlength was selected to be 60% that 
of the wing chordlength. Canard-wing horizontal 
and vertical separation distances were S = 1.0 
and D = 0.25, respectively, except where other- 
wise noted. 


Steady Transonic Interference 

Steady pressure distributions for the 
isolated airfoils are shown by the solid lines 
in Fig. 3. A shock wave of moderate strength Is 
present on the upper surfaces near 30% chord; 
the lower surfaces are entirely subcrltlcal. 
Steady pressure distributions on the closely- 
coupled canard-wing configuration are also shown 
in Fig. 3 (dashed lines). For this configura- 
tion, the canard produces a downwash on the wing 
which eliminates the shock such that the flow 
about the wing is entirely subcrltlcal. Con- 
versely, the wing produces an upwash on the 
canard such that the canard upper surface shock 
Is Increased In strength and located further aft 
near 44% chord. 

Transonic steady pressure distributions for 
a range of horizontal and vertical separation 
distances were obtained. Results for horizontal 
separation S = 0.85, 1.0, 1.25, 1.5, and 2.0, 
with constant vertical separation D = 0.25, are 
shown In Fig. 4. As the distance between canard 
and wing becomes small, the shock on the canard 
upper surface Increases in strength. The shock 
on the wing upper surface decreases In strength 
and then disappears, with decreasing horizontal 
separation distance. Results for a range of 
vertical separation distance show very similar 
steady transonic interference effects and are 
therefore not shown here. 



Mach number contour lines for the Isolated 
airfoils and several canard-wing configurations 
are shown In Fig. 5. The Mach contour Interval 
Is AM * 0.04. For the Isolated airfoils (Fig. 
5(a)), the shock wave Is represented by the 
close proximity of contour lines on the upper 
surfaces near 30% chord. For the canard-wing 
configurations (Figs. 5(b)-5(d)), the contour 
lines along the wing upper surface become 
relatively evenly spaced with the weakening and 
disappearance of the shock. The contour lines 
along the canard upper surface Indicate the 
presence of a strong shock wave especially 
for the closely-coupled configuration with S = 
1.0 and D = 0.25 (Fig. 5(d)). These Mach number 
contours clearly illustrate the steady transonic 
Interaction between the two airfoils. 


Unsteady Transonic Interference 

Unsteady transonic airloads were computed 
using the pulse transfer-function analysis 
method. Sample time histories for an Input wing 
pitch pulse and the resulting wing lift coeffi- 
cient and canard lift coefficient are shown In 
Fig. 6. Only the first 102 time steps are 
plotted. (Moment coefficient time histories are 
also computed for both canard and wing but 
are not shown here.) The aerodynamic forces In 



c; 

S=2.00 


S— 1 .50 


S«1.25 


S-1.00 


S-0.85 


Fig. 3 Steady pressure distributions on 
Isolated airfoils and canard-wing 
(S = 1.0, D = 0.25) NACA 0010 
configurations at M = 0.76 and 

a 0 = 1.0°. 


Fig. 4 Steady pressure distributions on NACA 
0010 canard-wing for a range of 
horizontal separation distance S 
(D » 0.25) at M = 0.76 and 

Oq = 1.0°. 
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the frequency domain are determined by dividing 
the FFT of the output lift coefficient time 
histories by the FFT of the input wing pitch 
pulse. The resulting wing lift coefficient due 
to wing pitch and canard lift coefficient due to 
wing pitch are shown In Figs. 7(a) and 7(b), 
respectively. These coefficients are plotted as 
real and imaginary functions of reduced fre- 
quency k. (The unsteady aerodynamic coeffi- 
cients Qij are defined as listed in Table 1.) 
To assess the applicability of the pulse 
transfer-function analysis to two-dimensional 






(c) S = 1.5, D = 0.25. 



(d) S = 1.0, D = 0.25. 


Fig. 5 Mach number contour lines (AM = 0.04) 
for isolated and various NACA 0010 
canard-wing configurations at M = 0.76 
and oq = 1.0°; 


Fig. 6 Wing pitch pulse and resulting wing 
and canard lift coefficient time 
histories for closely-coupled 
(S = 1.0, D = 0.25) NACA 0010 canard 
wing at M = 0.76 and a 0 = 1.0°. 




REDUCED FREQUENCY k REDUCED FREQUENCY k 


(a) wing lift coeffi- 
cient due to wing 
pitch. 


(b) canard lift coeffi- 
cient due to wing 
pitch. 


Fig. 7 Comparison between unsteady lift 
coefficients calculated by pulse 
analysis with oscillatory analysis for 
closely-coupled (S = 1.0, D = 0.25) 
NACA 0010 canard-wing at M =0.76 and 
a 0 = 1.0°. 
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Table 1.- Definition of unsteady aerodynamic 
coefficient Qjj. 


i 

j 

unsteady coefficient 

due to motion 

1 

1 

wing lift 

wing plunge 

1 

2 

wing lift 

wing pitch 

1 

3 

wing lift 

canard plunge 

1 

4 

wing lift 

canard pitch 

2 

1 

wing moment 

wing plunge 

2 

2 

wing moment 

wing pitch 

2 

3 

wing moment 

canard plunge 

2 

4 

wing moment 

canard pitch 

3 

1 

canard lift 

wing plunge 

3 

2 

canard lift 

wing pitch 

3 

3 

canard lift 

canard plunge 

3 

4 

canard lift 

canard pitch 

4 

1 

canard moment 

wing plunge 

4 

2 

canard moment 

wing pitch 

4 

3 

canard moment 

canard plunge 

4 

4 

canard moment 

canard pitch 




canard-wing configurations, harmonic oscillatory 
solutions were obtained for discrete values of 
reduced frequency. Unsteady wing and canard 
lift coefficients due to wing pitch, computed by 
the pulse transfer-function analysis, are com- 
pared with coefficients from the oscillatory 
analysis in Fig. 7. The close agreement between 
pulse and oscillatory forces thus verifies the 
extension of the pulse analysis to canard-wing 
geometries and demonstrates the ability of the 
pulse analysis to accurately predict transonic 
unsteady aerodynamic Interference effects. 

Self-induced unsteady airloads on the 
canard and wing due to motion of the respective 
airfoils, computed using the pulse tranfer- 
function analysis, are presented in Figs. 8 and 
9. These unsteady airloads are compared with 
the Isolated airfoil airloads to demonstrate the 
effects of aerodynamic Interference. Wing lift 
and moment coefficients due to wing plunge are 
shown in Figs. 8(a) and 8(b), respectively; wing 
lift and moment coefficients due to wing pitch 
are shown in Figs. 8(c) and 8(d), respectively. 
In general, the wing forces due to wing plunge 
motion show only small differences between 




(a) wing lift coeffi- 
cient due to wing 
plunge. 



(c) wing lift coeffi- 
cient due to wing 
pitch. 


(b) wing moment coeffi- 
cient due to wing 
plunge. 



(d) wing moment coeffi- 
cient due to wing 
pitch. 


(a) canard lift coeffi- 
cient due to canard 
plunge. 



REDUCED FREQUENCY k 


(c) canard lift coeffi- 
cient due to canard 
pitch. 


(b) canard moment coeffi- 
cient due to canard 
plunge. 



REDUCED FREQUENCY k 


(d) canard moment coeffi- 
cient due to canard 
pitch. 


Fig. 8 Unsteady wing coefficients for 
isolated airfoil and canard-wing 
(S = 1.0, D = 0.25) NACA 0010 config- 
urations at M = 0.76 and a 0 = 1.0°; 


Fig. 9 Unsteady canard coefficients for 
isolated airfoil and canard-wing 
(S = 1.0, D = 0.25) NACA 0010 config- 
urations at M = 0.76 and olq = 1.0°; 
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isolated airfoil and canard-wing cases. How- 
ever, large differences between the two sets of 
results occur in the real and imaginary parts of 
the wing lift coefficient due to wing pitch 
(Fig. 8(c)), for low values of reduced frequency 

k. Similar comparisons of self-induced unsteady 
airloads on the canard are shown in Fig. 9. 
Canard lift and moment coefficients due to 
canard plunge are shown in Figs. 9(a) and 9(b), 
respectively; canard lift and moment 
coefficients due to canard pitch are shown in 
Figs. 9(c) and 9(d), respectively. For the 
canard airloads due to canard motion, the 
differences between isolated airfoil and 
canard-wing results are generally small, except 
for the canard moment coefficient due to canard 
pitch (Fig. 9(d)). 

Oscillatory motion of the wing induces 
harmonic lift and moment on the fixed canard. 
Conversely, unsteady motion of the canard 
induces unsteady lift and moment on the fixed 
wing. Examples of these induced unsteady tran- 
sonic airloads as functions of reduced frequency 
are shown in Figs. 10 and 11. Results were 
obtained for a variety of horizontal and verti- 
cal separation distances, although unsteady 
forces are plotted here for S = 0.85, 1.0, 1.25, 

l. 5, and 2.0 with D = 0.25 only. Figures 10(a) 
and 10(b) show the canard lift coefficient due 
to wing plunge and the canard lift coefficient 
due to wing pitch, respectively. In both cases, 
the induced unsteady airloads become consider- 
ably larger when the separation distance becomes 
smaller as expected. Real and imaginary parts 
of the unsteady airloads, for a given separation 
distance S, are very similar in shape in compar- 
ison with those at other values of S, when the 
reduced frequency axis is either stretched or 
compressed. As shown in Fig. 10(a), for 
example, the first zero-crossing of the real 
part of Q31 occurs at successively increasing 
values of k for monotonical ly decreasing separa- 
tion distance S. The forces are also largest at 
low values of reduced frequency and tend to 
decrease in magnitude at higher values of 
reduced frequency. The large magnitudes at low 
reduced frequency are qualitatively consistent 
with the relatively long times associated with 
upstream propagating disturbances (from wing to 
canard) in comparison with downstream propagat- 
ing disturbances (from canard to wing). 

Figures 11(a) and 11(b) show the wing lift 
coefficient due to canard plunge and the wing 
lift coefficient due to canard pitch, respec- 
tively. In contrast with the results of Fig. 
10, these induced unsteady airloads do not tend 
to zero as the horizontal separation distance S 
becomes large. This is because the oscillatory 
wake of the canard is in close proximity (D = 
0.25) above the wing which influences wing pres- 
sures and hence wing lift, even for large values 
of S. Also in contrast with the results of 
Fig. 10, the induced airloads of Fig. 11 are of 
sizeable magnitude at the higher values of 
reduced frequency. This characteristic is 
attributed to the relatively shorter times for 
disturbances to propagate downstream from the 
oscillating canard to the motionless wing in 
comparison to upstream propagation from wing to 
canard. 



REDUCED FREQUENCY k 


(a) canard lift coefficient due to 
wing plunge. 



REDUCED FREQUENCY k 


(b) canard lift coefficient due to 
wing pitch. 


Fig. 10 Unsteady canard lift coefficients 

induced by wing motion for a range of 
horizontal separation distance S at 
M = 0.76 and a Q = 1.0° (D = 0.25); 
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Aeroelastic Applications 


r- CANARD 


WING 
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REDUCED FREQUENCY k 


(a) wing lift coefficient due to 
canard plunge. 



REDUCED FREQUENCY k 


(b) wing lift coefficient due to 
canard pitch. 


To investigate the effects of aerodynamic 
interference on aeroelastic stability and 
flutter. Pads model stability calculations were 
performed. Dynamic pressure root-loci of the 
wing for Isolated airfoil and closely-coupled 
canard-wing configurations are shown in 
Fig. 12. With increasing flight dynamic pres- 
sure q, the wing pitch mode moves to the, left in 
the stable left-half of the complex s-plane. 
The plunge dominated root-locus becomes the 
flutter mode at dynamic pressures of q £ = 
0.21 and q p = 0.35 for the isolated airfoil 
and closely-coupled canard-wing configurations, 
respectively. For the case considered, inclu- 
sion of the motionless canard in the transonic 
flowfield lowers the wing pitch modal frequency 
and increases damping in the wing plunge mode, 
thus delaying the onset of flutter. A 67% 
increase in flutter dynamic pressure q p (or 
equivalently, a 29% increase in flutter h speed 
index Vp) was obtained with the inclusion of 
the canard. The increase in wing flutter speed 
for the closely-coupled configuration is attri- 
buted to decreased transonic effects on the wing 
caused by the canard downwash. Flutter speeds 
for a range of canard-wing horizontal and verti- 
cal separation are plotted in Fig. 13. As the 
separation distance between the aeroelastic wing 
and the aerodynamical ly interfering canard 
decreases, the wing flutter speed increases. 
The flutter speed index Vp versus vertical 
separation distance D curve shown in Fig. 13 is 
not symmetric about the D = 0 line because the 
two airfoils are at 1.0° mean angle of attack. 
The flutter results would be symmetric about 
D * 0 if the airfoils were at 0° mean angle of 
attack. 
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Fig. 11 Unsteady wing lift coefficients 

induced by canard motion for a range 
of horizontal separation distance S 
at M = 0.76 and a 0 = 1.0° 

(D = 0.25); 


Fig. 12 Wing dynamic pressure root-loci 
for isolated airfoil and canard- 
wing (S = 1.0, D = 0.25) NACA 0010 
configurations at M = 0.76 and 
a 0 = 1.0°. 
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solving the three-dimensional canard-wing inter- 
action problem. These calculations were per- 
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sonic smal 1 -disturbance code to include an 
additional airfoil. The code is now capable of 
computing unsteady transonic flowfields about 
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geometries. 

For the configurations studied, the canard 
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weakened or eliminated the wing shock. 
Conversely, the wing produced an upwash field on 
the canard which increased the canard shock 
strength. Mach number contour lines for several 
canard-wing configurations clearly illustrated 
the steady transonic interaction between the two 
airfoils. Unsteady transonic airloads as a 
function of reduced frequency were computed by 
extending a pulse t ransfer-functi on analysis to 
canard-wing configurations. The accuracy of 
this analysis was confirmed by the excellent 
agreement found between unsteady airloads com- 
puted by pulse analysis and those calculated by 
oscillatory analysis. Results for a range of 
canard-wing horizontal separation distance 
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on transonic unsteady airloads. Aeroelastic 
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demonstrated the effects of aerodynamic inter- 
ference on aeroelastic stability and flutter. 
Inclusion of the canard in the transonic flow- 
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